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ABSTRACT 

Complex  temperature  anomalies  observed  near  Arctic 
ice  margins  in  the  Chukchi  Sea  were  found  to  be  associated 
with  the  interaction  of  the  warm  coastal  current  with  the 
ice  cover  and  resident  bottom  water.   These  anomalies  were 
characterized  by  large  temperature  gradients  found  in  areas 
of  very  low  density  gradient.   They  were  in  the  form  of 
inversions,  interleavings ,  and  discrete  parcels  which  varied 
greatly  in  short  distances.   The  processes  which  produced 
mesostructure  were  directly  linked  to  the  presence  of  ice 
and  were  found  to  affect  the  entire  water  column.   Structure 
was  found  in  the  vicinity  of  the  ice  margin  but  was  observed 
to  dissipate  well  inside  the  ice  margin.   A  dynamic  high  was 
found  which  was  related  to  the  melting  of  ice  and  character- 
istically occurred  in  the  melting  zone  of  the  ice  margin. 
This  dynamic  high  was  limited  to  a  depth  of  10  meters  and 
is  believed  to  be  a  factor  in  the  deepening  of  mesostructure 
elements  and  other  phenomena  as  they  pass  the  ice  margin 
with  the  current. 
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I.   INTRODUCTION 

A.  DESCRIPTION  OF  THE  PROBLEM 

An  investigation  was  conducted  to  examine  the  mechanisms 
associated  with  the  phenomenon  of  temperature  mesostructure 
found  in  the  water  columns  near  Arctic  ice  margins.   The 
term  mesostructure  describes  the  relatively  large  anomalies 
which  occur  in  the  vertical  temperature  profiles;  such  a 
profile  is  shown  in  Figure  1.   Corse  (1974)  defined  an 
element  of  mesostructure  as  that  portion  of  the  profile 
from  the  depth  where  the  temperature  gradient  became  positive 
to  the  depth  at  which  the  temperature  was  again  the  value 
it  had  when  the  gradient  became  positive.   These  anomalies 
were  considered  as  warm,  vice  cold,  since  they  were  found 
to  be  associated  with  the  advection  of  warm  water  into  the 
region.   The  magnitude  of  an  element  was  found  to  be  of 
the  order  of  meters  in  depth  and  of  tenths  of  degrees  in 
temperature  variation.   The  phenomenon  of  mesostructure 
was  found  to  be  associated  with  the  presence  of  an  ice 
margin  and  was  present  to  some  degree  in  all  ice  margin 
crossings  investigated  in  this  study. 

B.  DATA  BACKGROUND  AND  SOURCE 

Investigations  into  temperature  structure  near  ice 
margins  had  their  origin  in  reports  of  severe  deterioration 
of  sonar  propagation,  which  was  believed  to  be  the  result 
of  complex  sound-speed  profiles  near  the  ice  margins.   As 
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a  result  of  this  interest,  research  cruises,  under  the 
general  project  title  MIZPAC,  were  directed  by  the  Arctic 
Submarine  Laboratory,  Naval  Undersea  Center.   These  cruises, 
carried  out  in  1971,  1972,  and  1974  on  Coast  Guard  ice 
breakers,  have  resulted  in  a  data  base  of  350  STD  stations 
with  another  such  cruise  scheduled  for  1975.   The  results 
of  the  1971  and  1972  cruises  were  published  (Paquette  and 
Bourke,  1973,  1974)  and  those  for  1974  and  1975  are  in 
preparation  by  Paquette  and  Bourke.   The  data  from  the  1974 
cruise  were  analyzed  because  they  were  the  most  useful  of 
the  available  data,  due  to  the  number  of  ice  crossings  and 
the  close  proximity  of  the  stations. 

The  1974  data  were  collected  using  a  Bissett-Berman 
Model  9006  STD.   Two  lowerings  were  required  for  each 
station;  a  deep  cast,  and  a  shallow  cast  with  the  conductivity 
sensor  shunted  with  a  precision  resistor  to  enable  recording 
of  the  low  salinities  encountered  near  the  surface.   Seven 
series  of  closely  spaced  stations  were  analyzed,  each  one 
representing  a  penetration  of  the  ice  margin  and  termed  a 
"crossing."  Figure  2  indicates  the  location  of  the  crossings, 
carried  out  during  the  period  18  July  to  26  July  1974,  and 
several  positions  of  the  ice  margin  for  this  period.   The 
positions  of  the  ice  margin  were  derived  from  cruise  data 
and  NOAA-4  VHRR  satellite  photographs  interpreted  by  the 
NOAA-NESS  Environmental  Products  Group. 
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C.   DESCRIPTION  OF  THE  MIZPAC  74  AREA 

The  Chukchi  Sea  is  a  marginal  sea  with  very  little 
bottom  relief  in  the  area  of  interest  and  with  an  average 
depth  of  about  45  m. 

The  ice  coverage  and  thus  the  ice  margin  varied 
considerably  during  the  months  of  June  through  August 
1974  (Fig.  2).   There  is  much  year-to-year  variation.   In 
1974  the  ice  melted  back  from  Bering  Strait  in  June  to 
73°N  by  the  end  of  August,  farther  north  than  normal. 

The  pattern  of  the  melting  process  followed  closely 
that  of  the  previously  known  current  system  of  the  Chukchi 
Sea.   The  summer  surface  circulation  of  the  Chukchi  Sea  is 
dominated  by  a  general  northward  flow  of  warm  water  from 
Bering  Strait.   This  flow  was  observed  by  Paquette  and 
Bourke  (1974)  and  was  typified  by  a  surface  temperature  of 
5°  to  11°C  during  late  June  through  September  and  probably 
involved  the  upper  mixed  layer  of  10  to  20  m  depth.   Little 
is  known  about  the  movement  of  the  lower  layer.   The  current 
was  described  as  staying  to  the  right  side  of  the  Chukchi 
Sea  and  splitting  at  Cape  Lisburne,  one  branch  following 
the  Alaskan  coast  and  the  other  turning  north-north-west 
toward  Herald  Shoal.   Garrison  et  a.  (1973)  reported  the 
speed,  relative  to  an  ice  floe,  of  the  current  along  the 
Alaskan  coast  near  72°N  to  vary  between  10  and  50  cm  sec   , 
with  a  velocity  shear  between  10  and  20  m  which  was  usually 
associated  with  the  pynocline. 
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D.   OBJECT  OF  THE  STUDY 

The  thrust  of  this  undertaking  was  to  gain  an  under- 
standing of  the  mechanisms  by  which  mesostructure  is  formed, 
utilizing  new  data  and  previous  findings.   It  was  particu- 
larly directed  toward  finding  the  relationship  between  the 
formation  of  mesostructure  and  the  following  processes: 
the  influx  of  warm  water  and  the  melting  of  the  ice  cover, 
the  sloping  of  the  density  surfaces,  the  dynamic  heights 
near  the  ice  margin,  and  local  currents.   Additionally, 
it  was  desired  to  determine  if  different  processes  were 
responsible  for  deep  and  shallow  mesostructure. 
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II.   INVESTIGATIONAL  PROCEDURE 

The  data  presentations  were  predicated  on  the  assumption 
that  warm  water  flowed  northward  toward  the  ice  and  that  the 
interaction  of  ice  and  water  resulted  in  mesostructure .   It 
was  expected  that  the  phenomena  involved  would  be  evident 
in  vertical  sections  of  the  conventional  oceanographic 
properties  along  the  lines  of  the  several  ice-margin 
crossings.   Temperature  and  sigma-t  sections  were  therefore 
prepared.   The  latter,  being  closely  related  to  salinity 
at  locally  observed  water  temperatures,  eliminated  the 
need  for  salinity  sections. 

Later  it  became  apparent  that  the  question  of  direction 
of  flow  of  the  warm  water  could  be  investigated  by  means  of 
the  heat  and  ice  melt  contents  of  the  water  column.   The 
local  dynamic  height  gradient  near  the  ice  margin  seemed 
a  possible  driving  force  for  downward  mixing.   These 
parameters  also  were  computed  for  the  ice-margin  crossings. 
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III.   DISCUSSION 

A.   TEMPERATURE  CHARACTERISTICS 

The  temperature  profiles,  nested  in  sequence,  for  each 
ice  margin  crossing  are  shown  in  Figures  3  through  9.   The 
station  number  and  bottom  temperature  (°C)  are  indicated 
at  the  bottom  of  each  trace,  and  the  ice  concentration  in 
oktas  (eighths)  is  displayed  at  the  top  with  the  symbol 
MX"  indicating  no  data.   The  station  numbers  increased  as 
the  ice  margin  was  penetrated  from  outside  the  ice;  thus 
the  crossings  proceeded  from  left  to  right  in  all  figures. 

The  temperature  profiles  displayed  several  characteristics 
which  were  common  to  all  the  crossings  in  general.   A  layer 
of  relatively  warm  water,  between  the  surface  and  20  m  depth, 
appeared  outside  the  ice  in  all  cases.   This  layer  was 
cooled  and  modified  with  progress  toward  and  across  the  ice 
margin  and  was  greatly  diminished  or  absent  well  inside  the 
ice  margin.   Surficial  cooling  caused  the  temperature 
maximum  of  this  layer  to  descend  several  meters  when  ice 
was  encountered  forming  the  "nose"  feature  described  by 
Paquette  and  Bourke  (1973) . 

Corse  (1974)  described  mesostructure  as  being  divided 
into  three  classes:   the  nose  feature,  and  shallow  and  deep 
structure.   He  defined  shallow  structure  as  having  sigma-t 
values  less  than  25.5  and  deep  structure  as  that  with 
greater  values.   He  placed  the  sigma-t  value  of  25.5  at 
about  20  m  depth;  however,  as  will  be  seen,  this 
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investigation  found  the  26.0  surface  at  about  20  m.   This 
depth  was  used  to  make  a  distinction  between  shallow  and 
deep  mesostructure . 

Shallow  structure  was  present  in  varying  degrees  of 
severity  in  all  crossings.   Crossing  3  exhibited  the  most 
severe  shallow  structure  as  compared  to  the  lesser  intensi- 
ties found  in  Crossings  4,  5,  and  6  and  to  Crossings  1  and 
2,  which  displayed  very  little  shallow  structure.   The 
presence  of  shallow  structure  correlated  well  with  the 
presence  of  the  nose  feature,  which  was  present  in  all 
crossings  except  1  and  2.   Thus,  the  nose  feature  appeared 
to  have  a  causal  relationship  with  shallow  structure. 

Deep  structure  was  found  to  be  severe  in  Crossings  1 
and  2  and  very  slight  or  absent  in  the  other  crossings. 
The  presence  of  deep  structure  appeared  to  be  a  separate 
phenomenon  from  that  of  shallow  structure. 

The  temperature  profiles  displayed  no  apparent 
correlation  between  the  presence  of  mesostructure  and  the 
maximum  temperature  present.   This  differed  from  the 
belief  of  Corse  (1974)  in  the  relationship  between  high 
temperatures  and  the  presence  of  deep  mesostructure. 
However,  the  formation  and  presence  of  structure,  both 
shallow  and  deep,  did  correlate  with  the  presence  of  an 
ice  margin,  as  did  the  formation  and  descent  of  the  nose 
feature.   The  mesostructure  eventually  almost  disappeared 
at  some  point  within  the  ice  margin.   This  may  indicate 
that  the  structure  had  been  dissipated  as  it  moved 
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"downstream"  with  a  northerly  flowing  current  past  the  ice 
edge  or  that  there  was  little  or  no  water  flow  toward  the 
ice. 

Vertical  cross-sections  of  temperature  for  the  seven 
crossings  are  presented  in  Figures  10  through  16.   Again 
the  relatively  warm  surface  layer,  with  a  maximum  tempera- 
ture of  3  to  6°C,  may  be  observed  south  of  the  ice  at 
depths  of  less  than  20  m.   This  layer  cooled  drastically  and 
lost  most  of  its  warmth  in  the  vicinity  of  the  ice  margin. 
This  was  obviously  due  to  a  utilization  of  this  heat  for 
the  melting  process  which  will  be  discussed  in  a  later 
section.   Referring  to  the  cross-sections,  the  isotherms 
were  very  irregular  and  displayed  large  gradients,  many 
inversions,  interleavings ,  and  lenses.   These  irregularities 
were  more  pronounced  in  the  vicinity  of  the  ice  margin  and 
dissipated  when  well  inside  it,  as  has  already  been  stated 
in  the  description  of  the  profiles.   Individual  elements  of 
mesostructure  usually  could  not  be  traced  for  more  than  two 
consecutive  stations,  a  distance  of  1  to  5  km.   The  maximum 
distance  over  which  a  mesostructure  feature  was  traced 
was  about  11  km.   This  is  in  reasonable  agreement  with  the 
results  of  Corse  (1974)  .   He  found  typical  correlation 
distances  of  1  to  2  km.   Structure  was  found  to  exist  for 
some  time  after  the  disappearance  of  ice  from  an  area 
since  it  was  found  as  far  as  28  km  south  of  the  ice  margin 
where  the  ice  had  apparently  been  recently. 
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Also  well  correlated  with  the  ice  margin  was  the  presence 
of  cold  bottom  water  which  was  absent  to  the  south  of  the 
margin.   This  type  of  water  appeared  fairly  precisely  at  the 
margin  and  thickened  as  one  progressed  into  the  ice.   It 
had  a  temperature  of  -1.6  to  -1.7°C  and  was  usually  either 
in  the  form  of  a  wedge,  as  in  Crossings  1,  2,  3,  and  7,  or 
that  of  a  cold  tongue  at  about  30  m  depth  as  in  Crossings  4, 
5,  and  6.   This  correlation  of  the  boundary  of  the  cold 
bottom  water  with  the  ice  margin  is  indicative  of  the  trans- 
fer of  heat  downward  in  the  water  column  near  the  margin. 

The  dissipation  of  the  surface  heat  across  the  margins, 
the  melting  of  the  ice,  and  the  absence  of  the  cold  bottom 
water  south  of  the  margin  all  indicate  large  heat  transfers 
throughout  the  entire  water  column  in  the  vicinity  of  the 
ice  margin.   The  extreme  temperature  gradients  associated 
with  the  presence  of  structure  were  found  predominantly 
near  the  ice  margin,  but  not  far  under  the  ice,  indicative 
of  strong  downward  mixing  processes  in  the  vicinity  of  the 
margin.   Farther  back  under  the  ice,  however,  this  powerful 
downward  mixing  of  near-surface  water  appeared  to  be  lacking. 
Therefore,  ice  keels  probably  are  not  the  cause  of  mixing 
as  was  suggested  by  Corse  (1974) .   Instead,  the  mixing  must 
be  due  to  a  dynamic  process  which  is  peculiar  to  the  ice 
margin.   The  process  must  be  due  to  something  distinctive 
of  the  ice  margin,  most  likely  the  rapid  melting  of  ice. 
This  melting  process  and  its  relationship  to  the  formation 
of  mesostructure  will  be  discussed  in  a  subsequent  section. 
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B.   DENSITY  CHARACTERISTICS 

The  cross-sectional  density  distributions,  plotted  in 
sigma- t  units,  are  shown  in  Figures  17  through  23.   All 
crossings  were  typified  by  a  high  gradient  layer  in  the 
upper  10  m,  with  localized  pockets  of  low  density  water 
located  at  the  surface  which  correspond  to  areas  of  melting 
ice.   The  density  surfaces  like  the  isotherms,  displayed 
irregularity  in  the  vicinity  of  the  ice  margin.   However, 
the  density  surfaces,  unlike  the  isotherms,  displayed  no 
inversions.   The  density  at  the  surface  was  observed  to  be 
higher  outside  the  margin  and  decreased  near  it.   These 
phenomena  were  undoubtedly  due  to  melting  processes.   The 
irregularity,  high  density  gradients,  and  low  density 
pockets  did  not  appear  farther  back  under  the  ice.   The 
rapidly  added  melt  water  must  have  been  mixed  downward 
rather  quickly  in  the  vicinity  of  the  margin  so  as  to  have 
decreased  the  high  gradients  that  were  present.   The  density 
cross  sections  display  the  same  general  characteristics  as 
do  the  temperature  cross  sections  and  there  is  evidence  of 
dynamic  activity  in  all  cases. 

The  deepest  isopycnals,  present  well  inside  the  ice 
margin,  were  destroyed  by  downward  mixing  of  melt  water 
near  the  margin.   These  isopycnals  usually  rose  toward  the 
interior  of  the  ice  field,  forming  a  wedge-like  feature. 
In  Crossings  1,  2,  3,  and  4  the  wedge  was  pronounced.   In 
Crossings  5  and  6  the  wedge  was  less  pronounced  and  did 
not  exist  in  Crossing  7.   This  wedge  was  well  correlated 
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with  the  cold,  bottom  water  found  in  the  temperature 
data. 

The  complex  tilting  of  the  density  surfaces  is  evidence 
of  lateral  pressures  that  must  be  maintained  by  dynamic 
processes.   Thus,  stronger  density  structures  should  be 
indicative  of  strong  dynamic  processes  and  more  severe 
mesostructure .   Deep  mesostructure  was  found  to  be  more 
severe  for  those  crossings  with  a  well-defined  wedge  and 
thus  gave  support  to  this  hypothesized  relationship  between 
deep  mesostructure  and  strong  dynamic  processes.   This 
relationship  was  further  confirmed  by  an  investigation  of 
the  dynamic  heights  associated  with  the  marginal  ice  zone 
and  is  discussed  in  a  subsequent  section. 

Observations  of  the  26.0  and  26.5  sigma-t  surfaces 
revealed  that  they  had  a  pattern  of  variable  spacing 
related  to  the  ice  edge.   For  Crossings  1,  2,  3,  and  4 
(Figures  17-20)  these  isopycnals  were  more  widely  separated 
outside  the  ice  margin,  thus  defining  a  zone  of  very  low 
vertical  density  gradient.   For  Crossings  5,  6,  and  7 
(Figures  21-23)  this  low  density  gradient  zone  was  found 
well  inside  the  ice  edge,  but  in  Crossing  5  it  occurred 
between  the  26.5  and  26.8  sigma-t  surfaces.   Inspection  of 
the  above  figures  indicates  that  most  of  the  deep  meso- 
structure is  found  within  these  zones  of  weak  density 
gradient.   This  .is  perhaps  an  expected  phenomenon  as 
relatively  little  energy  is  required  to  mix  anomalous 
water  downward  in  the  water  column  when  the  density 
gradient  is  small. 
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It  is  important  to  remember  at  this  point  that  a 
temperature  or  density  cross  section  suggests  that  there 
is  a  transport  of  properties  from  left  to  right.   This  need 
not  be  the  case;  flow  may  be  oblique  or  nearly  parallel  to 
the  ice  edge.   This  type  of  flow  pattern  may  account  for 
small  degrees  of  tilt  to  the  isopycnals  within  the  ice 
or  for  weak  dynamics  and  weak  structure  to  exist.   It  may 
also  cause  a  low  rate  of  heat  transport  to  the  ice  and 
consequently  low  melting  rates.   Considering  these  factors 
the  heat  balance  across  the  ice  margin  was  investigated. 

C.   HEAT  BALANCE  ACROSS  AN  ICE  MARGIN 

The  disappearance  of  the  warm  surface  water  under  the 
ice  margin  and  the  melting  which  occurred  led  to  an  analy- 
sis of  the  heat  content  of  the  water  columns  associated  with 
the  ice  margin  crossings.   Specifically,  it  was  desired  to 

correlate  the  amount  of  heat  dissipated  at  the  margin  with 

< 

the  melting  of  an  equivalent  amount  of  ice  and  to  test  the 
associated  idea  that  water  flows  along  the  section  toward 
the  ice.   If  it  did  not,  it  might  be  an  explanation  for 
weak  dynamic  mixing  processes  and  weak  structure. 

The  heat  content  of  the  water  column  at  each  station 
in  the  various  crossings  was  calculated,  using  as  a  ref- 
erence water  of  33.33°/oo  and  -1.8°C.   This  reference  was 
selected  since  33.33°/0o  is  the  salinity  of  freezing 
sea  water  at  the  lowest  water  temperature  observed  in  the 
area.   The  resultant  values  were  termed  the  "heat  excess" 
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as  defined  by  equation  (1) 


-D 


H  =  C 


p(T-TR)  dz 


(1) 


where:    H  =  Heat  excess  (cal-cm   ) 

C   =  Specific  heat  (cal-gm   °C   ) 

_  3 
p  =  Density  (gm-cm   ) 

T  =  Temperature  (°C)  at  depth  z  cm 

TR  =  Reference  temperature  (°C) 

D  =  Depth  of  the  water  column  (cm) 

dz  =  Depth  increment  (cm) 

In  order  to  examine  the  melting  phenomenon  in  terms  of 

the  fresh  water  introduced  from  melting  ice  an  "ice  melt 

content"  was  also  defined.   It  was  defined  as  that  weight 

of  ice,  with  an  initial  salinity  of  6.0°/oo  and  a  density 

-3 
of  0.94  gm  cm   ,  that  must  be  transformed  into  melt  water 

to  produce  the  salinity  structure  present  in  the  given 

column  of  water,  referenced  again  to  a  water  column  having 

an  initial  salinity  of  33.33°/00.   The  ice  melt  content  is 

defined  by  equation  (2)  : 


where : 


I  = 


rB  sR  -  s 

p(5T— 5")  dz 


J 


0 


R 


I 


I  = 

SR  " 
S  = 

ST  = 


-2. 


Ice  melt  content  (gm-cm   ) 
Reference  salinity  (ppt) 
Salinity  (ppt)  at  depth  z  cm 
Salinity  of  ice  (ppt) 


(2) 
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Values  of  0.956  cal  gm"1  0C_1  for  C   and  100  cm  for  dz 

P 

were  selected  and  the  values  for  the  remaining  variables 
(T,  D,  p,  and  S)  were  obtained  from  the  individual  station 
data.   Using  these  values  the  total  integrated  heat  excess 
and  ice  melt  content  from  the  surface  to  the  bottom  were 
computed.   The  values  from  the  surface  to  35  m  were  selected 
for  comparison  purposes  for  the  analysis.   One  may  make 
nearly  as  good  an  argument  for  integrating  to  the  bottom. 
However,  the  likelihood  of  incorporating  an  anomalous 
contribution  from  a  variable  depth  limit  led  to  the  choice 
of  a  fixed  limit.   The  results  of  these  computations  for 
all  crossings  are  displayed  in  Table  I. 

Ideally,  with  a  flow  of  water  directly  across  the  ice 
margin,  the  heat  excess  and  ice  melt  content  data  should 
present  an  inverse  relationship.   The  heat  excess  should 
decrease  and  remain  low  after  passing  the  margin  as  the 
ice  melt  increases,  which  in  turn  should  remain  high  due  to 
the  melting  process.   Considering  the  heat  excess  and  ice 
melt  content  data,  many  irregularities  were  present  in  each 
of  the  crossings.   However,  all  crossings  did  display  a 
general  heat  decrease  and  an  increase  in  ice  melt  content 
close  to  the  ice  margin,  except  for  Crossing  6  which 
displayed  a  decrease  in  ice  melt  content  and  heat  excess. 
Crossing  2  is  presented  as  an  example  (Figure  24) .   It 
was  close  to  the  average  of  all  the  crossings,  excluding 
Crossing  6,  even  to  the  extent  of  being  less  irregular  than 
most.   The  data  from  the  seven  crossings  revealed  a  decrease 
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-2 
in  excess  heat  ranging  from  4900-10300  cal  cm   across 

_  2 
the  ice  margin.   The  average  change  was  6600  cal  cm 

Equation  (3a)  describes  the  heat  (AH)  that  is  required 
to  warm  and  melt  ice  of  initial  temperature  T  and  thick- 
ness h 


AH  =  [(T   -  T  )  C   +  AT]hp  (3a) 


_2 
where:     AH  -  Heat  change  (cam  cm  ) 

T  =  Initial  temperature  (°C) 

T  =  Melting  temperature  (-.053S)  (°C) 

S  =  Salinity  (ppt) 

C  =  Specific  heat  of  fresh  ice  (cal  gm   °C   ) 

A-,  =  Latent  heat  of  fusion  of  sea  ice  (cal  gm  ) 

h  =  Ice  thickness  (cm) 

Solving  for  the  ice  thickness,  h,  we  arrive  at  equation  (3b) 


h  =  —  [(T   -  T  )  C   +  AT]_1.  (3b) 

P   L  v  g    oJ      p    IJ  K      J 


An  ice  cover  of  salinity  6.0°/oo,  initially  at  -2°C, 
and  with  a  density  of  0.94  gm  cm   was  selected  as  repre- 
sentative of  the  area  of  study.   This  selection  resulted 

in  values  of  0.5  cal  gm   °C   for  C   and  67  cal  gm   for 

-2 
X-r.      Using  the  average  change  of  heat  content,  6600  cal  cm 

across  the  margin  and  the  above  values,  an  ice  thickness 

of  1.0  m  would  be  melted. 
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The  ice  melt  content,  with  one  exception,  displayed  an 

-  2 
increase  across  the  margins  of  30-73  gm  cm   with  an 

-2 
average  value  of  54  gm  cm   .   The  exception  was  Crossing  6 

which  resulted  in  a  negative  value.   The  average  value  of 

-2 
54  gm  cm   when  equated  to  a  thickness  of  reference  ice  by 

equation  (4)  resulted  in  an  equivalent  ice  thickness  melted 

of  0.6  m: 

h  =  Alp'1  (4) 

where:    h  =  Thickness  of  ice  melted  (cm) 

_  2 
AI  =  Change  in  ice  excess  (gm-cm   ) 

_  3 
p  =  Density  of  ice  (gm-cm   ) . 

The  values  of  the  thickness  of  ice  melted  computed 
from  the  heat  excess  and  ice  melt  content  agree  favorably 
with  the  ice  thickness  of  about  0.8  m  predicted  for  the 
region  of  study  by  Wittman  et  al .  (1963).   This  lends 
support  to  the  assumptions  that  the  dissipation  of  the  heat 
across  the  ice  margin  can  be  accounted  for  by  the  melting 
of  the  ice  cover  and  that  the  flow  under  the  ice  margin 
is  generally  northward  (left  to  right  in  the  sections) . 

The  overall  changes  in  heat  excess  and  ice  melt  content 
across  the  margin  were  not  apparently  related  to  the  severity 
of  structure  present  (Table  III).   However,  the  rates  at 
which  they  changed  near  the  margin  were  possibly  related 
to  the  severity  of  shallow  mesostructure .   Crossings  1 
and  2  had  the  least  rapid  heat  and  melt  water  changes  in 
the  vicinity  of  the  margin  and  also  displayed  the  least 
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severe  shallow  structure.   Crossings  3,  4,  and  5  had 
significantly  more  rapid  heat  exchanges  at  the  margin  and 
also  had  the  more  severe  shallow  structure.   These  factors 
suggest  a  correlation  between  the  rate  of  heat  exchange 
and  the  severity  of  shallow  structure. 

In  most  cases  the  major  portion  of  the  heat  exchange 
occurred  at  the  ice  margin;  however,  in  Crossings  1  and 
2  it  occurred  outside  the  ice  margin.   This  can  be  explained 
by  a  flow  away  from  the  ice  subsequent  to  a  period  of 
melting  or  by  a  retreat  of  the  ice  front  that  was  more 
rapid  than  the  northward  flow  of  water. 

In  Crossings  1,  3,  4,  5,  and  7  there  was  an  irregular 
increase  of  ice  melt  content  toward  the  north  which  is 
what  one  might  expect  if  the  melt  water  near  the  margin 
were  carried  under  the  ice  by  the  current.   Crossings  2 
and  6  displayed  a  gradual  decrease  of  ice  melt  content 
as  one  progressed  farther  under  the  ice.   Possible  explana- 
tions for  this  difference  are  that  the  axis  of  flow  of  the 
northward  current  was  not  directly  across  the  margin  or 
an  irregularity  of  flow  occurred.   Thus,  it  appears  that 
the  flow  at  the  margins  is1 not  as  simple  as  was  assumed, 
although  it  is  generally  northward. 

D.   DYNAMIC  HEIGHTS  AND  LOCAL  CURRENTS 

The  degree  of  disturbance  in  the  density  surface 
(Figures  17-23)  and  the  presence  of  many  factors  which 
pointed  to  the  existence  of  highly  dynamic  processes  acting 
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in  the  immediate  vicinity  of  the  ice  edge  gave  rise  to  an 
interest  in  the  dynamic  heights  associated  with  melting  at 
the  marginal  ice  zone.   Neshyba's  model  (1974)  of  the  pro- 
grading  ice  cover  during  the  freeze-up  of  the  Bering  Sea 
predicted  the  inverse  phenomenon,  a  dynamic  low  under  the 
ice  margin.   His  model  described  the  freezing  process 
which  rejected  brine  and  mixed  it  downward  through  the  entire 
water  column  resulting  in  higher  density,  a  dynamic  low,  and 
a  horizontal  pressure  gradient  toward  the  ice  under  the 
freezing  zone.   Neshyba  concluded  that  freezing  rates  were 
slow  enough  such  that  the  pressure  gradient  must  be  in  geo- 
strophic  balance  and  must  therefore  correspond  to  a 
geostrophic  flow  to  the  east  under  the  freezing  zone.   His 
assumption  that  the  oceanic  system  was  otherwise  static 
probably  led  to  unrealistically  high  values  of  salinity 
and  correspondingly  high  negative  values  of  dynamic  height. 
An  extension  of  Neshyba's  work  to  the  melting  season  would 
suggest  that  the  converse  situation,  a  dynamic  high,  might 
be  found  under  a  melting  ice  zone  and,  with  melting  much 
more  rapid  than  freezing,  that  the  change  in  geopotential 
should  be  much  greater. 

Dynamic  heights  for  all  crossings  were  computed  utilizing 
the  HYDRO  program  (W.  R.  Church  Computer  Center,  NPS)  and 
the  results  appear  in  Table  II.   Figure  25  presents  a 
vertical  cross-section  of  dynamic  heights  for  Crossing  2 
which  was  selected  to  illustrate  significant  features.   All 
crossings  in  general  displayed  local  irregularities  and 
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all,  with  the  exception  of  Crossing  6,  displayed  a  hill 
(or  half  of  a  hill  as  in  Crossings  1,  3,  and  4)  situated 
near  the  ice  margin.   Evidence  of  this  dynamic  hill  was 
observed  to  a  depth  of  about  10  m  in  Crossings  1,  2,  3,  and 
4  and  to  about  15  m  in  Crossings  5  and  7.   Thus  the  hill 
was  generally  noticeable  in  the  near  surface  zone  only. 

These  results  were  somewhat  as  predicted.   However, 
the  hill  was  found  to  be  a  localized  phenomenon  extending 
over  horizontal  distances  of  up  to  30  km  across  the  margin 
and  its  effects  were  observed  to  only  about  10  m  of  depth. 
This  was  different  from  Neshyba's  freezing  model  which 
resulted  in  a  general  depression  associated  with  the  entire 
freezing  zone  which  extended  100  km  back  under  the  ice, 
and  whose  effects  were  notable  throughout  the  entire  water 
column.   The  heights  of  these  hills  were  of  the  order  of 
2  dyn  cm  as  compared  to  depressions  of  2.7  dyn  cm  calcu- 
lated by  Neshyba.   The  hills  due  to  melting  were  expected 
to  have  been  several  times  as  large  as  the  depressions  due 
to  freezing,  but  were  not.   This  was  probably  because  the 
water  at  the  ice  margin  was  not  static  as  in  Neshyba's 
model.   It  will  be  seen,  however,  that  the  magnitude  of 
this  hill  is  not  the  critical  factor. 

The  author's  data  revealed  that  the  height  of  the  hill 
under  the  ice  apparently  was  not  correlated  with  the  severity 
of  the  mesostructure .   However,  the  magnitude  of  the  hori- 
zontal gradient  of  the  dynamic  height  at  the  margin  did 
appear  to  be  related  directly  to  the  severity  and  presence 
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of  deep  mesostructure .   Table  III  shows  that  Crossings  1, 
2,  3,  and  7  had  the  highest  dynamic  height  gradients  and 
also  the  most  severe  deep  structure,  whereas  Crossings  4 
and  5  had  relatively  small  gradients  and  either  lacked 
or  had  very  little  deep  structure. 

The  gradient  observed  in  the  dynamic  heights  must 
result  in  a  lateral  pressure  gradient.   This  gradient  is 
probably  of  recent  origin  due  to  the  rapid  retreat  of  the 
ice  margin;  thus  geostrophic  balance  is  probably  not  achieved 
and  the  pressure  gradient  is  most  likely  balanced  against 
inertial  acceleration  or  internal  friction.   In  either  case 
the  pressure  gradient  leads  to  a  velocity  shear  in  the 
vertical  plane  normal  to  the  ice  margin.   Because  the 
pressure  gradients  on  the  two  opposing  sides  of  the  dynamic 
hill  are  oppositely  directed,  there  will  be  a  diminution 
of  flow  near  the  surface  on  the  seaward  side  of  the  dynamic 
hill  and  an  acceleration  of  flow  on  the  other.   Continuity 
of  flow  then  would  demand  that  the  shear  at  the  surface 
be  compensated  by  a  reverse  shear  at  deeper  depths  under 
the  dynamic  high.   It  is  believed  that  this  downward  directed 
flow  may  be  the  agency  which  caused  the  mixing  at  depth 
which  was  observed  near  the  ice  margin. 

Both  of  the  possible  opposing  forces  of  friction  and 
inertia  were  found  to  be  of  the  same  order  of  magnitude, 
depending  upon  one's  choice  of  viscosity  and  characteristic 
velocity,  but  viscosity  is  likely  to  exert  the  largest 
effect.   Therefore,  to  explore  the  magnitudes  of  the 
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effects,  approximate  calculations  were  made  based  on  the 
assumption  that  the  lateral  pressure  gradient  was  balanced 
by  the  frictional  shear  of  eddy  viscosity.   A  model  was 
developed  assuming  a  current,  u,  outside  the  ice,  flowing 
in  the  plane  of  the  pressure  gradient  and  uniform  in 
velocity  throughout  the  upper  warm  layer. 

The  step-like  profile  of  Figure  26  (solid  line)  was 
designed  for  the  current  velocity  structure  hypothesized 
outside  the  ice  margin  based  on  observed  currents  from 
other  investigations.   Equation  (5)  was  used  to  describe 
the  balance  of  frictional  shear  with  the  horizontal  pressure 
gradient : 

A  SjJ.    AD  ro 

p  3z2  "  AL  ^ 

A  =  Coefficient  of  eddy  voscosity  (gm-cm   sec   ) 

_  3 
p  =  Density  (gm-cm  ) 

u  =  Current  speed  (cm-sec   ) 

z  =  Depth  (cm) 

jj-   =  Horizontal  gradient  of  dynamic  height  (dyn  m-km   ) 

Equation  (5)  was  solved  for  u  with  the  conditions  that  the 

dynamic  effect  was  limited  to  10  m  and  that  z  =  10  m, 

—  =    0  and  u  =  15  cm  sec"  .   Specification  of  this  velocity 

o  Z 


was  not  necessary  to  the  solution,  but  it  did  influence  to 

s 

which  required  that  p  and  A  also  be  constants. 


3  2u 
ome  extent  the  choice  of  A.   £ — T   was  assumed  to  be  constant 

9  z 
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Figure    26.      Effect   of   the    Dynamic   Hill 
on   the    Current    Profile. 
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The  values  of  the  horizontal  gradient  of  dynamic 
height,  which  were  observed  for  the  crossings,  ranged 
between  0.0010  and  0.0039  dyn  m  km   with  an  average  of 
0.0026  dyn  m  km   for  all  crossings  except  Crossing  6 
(Table  III).   Values  for  density  and  depth  were  selected  as 
1.02  gm  cm   and  10  m,  respectively.   Selecting  a  value  for 
the  coefficient  of  eddy  viscosity  (A)  was  more  difficult. 
Eddy  viscosity  coefficients  have  an  extreme  range  of  values 
which  vary  with  the  scale  and  rate  of  the  phenomenon  in 

question.   Neuman  (1966)  shows  values  of  from  1  to  3000 

-1     - 1 

gm  cm   sec   .   Campbell  (1965) ,  however,  arrived  at  a 

-1     -1 
value  of  55  gm  cm   sec   for  a  numerical  model  of  the  ice 

packs  of  the  Arctic  Ocean,  in  which  a  relative  velocity  of 

13  cm  sec   existed  between  an  ice  floe  and  the  sea  water. 

Figure  27  shows  how  a  variation  in  A  affects  the  magnitude 

of  the  total  velocity  diminution  experienced  at  the  surface 

due  to  the  shear  forces  opposed  to  the  flow.   Varying  A 

from  50  to  500  gm  cm"   sec   reduces  the  current  from  23.5 

to  2.35  cm  sec   which  is  of  the  proper  order  of  magnitude 

for  the  currents  assumed  for  the  region.   An  intermediate 

-1     -1 
value  of  100  gm  cm   sec   resulted  in  a  decrease  of  the 

velocity  at  the  surface  of  13  cm  sec   ,  a  fortuitous  agree- 
ment with  Campbell. 

The  application  of  this  velocity  shear  in  the  upper  10  m 
of  the  water  column  and  the  enforcement  of  continuity  of 
flow  upon  the  initial  current  profile  of  Figure  26  resulted 
in  the  second  profile  (dashed  line)  of  Figure  26,  the 
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COEFFICIENT    OF   EDDY    VISCOSITY      (gmcm'Sec1) 


Figure  27.   Effect  of  Eddy  Viscosity  on 
Current  Velocity  Diminution. 
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profile  for  the  flow  beneath  the  ice  margin.   This  profile 
displays  a  descent  of  the  core  of  the  current  as  it 
encounters  the  margin  and  its  effects.   The  compensating 
flow  at  depth  is  pictorial;  it  did  not  seem  feasible  to 
calculate  it  in  a  more  meaningful  way.   This  restriction  to 
the  flow,  the  resultant  downward  motion,  and  the  accelera- 
tion of  flow  at  greater  depths  could  provide  a  mechanism 
for  the  mixing  process  which  occurred. 

It  is  evident  that  where  dynamic  height  gradients  are 
large,  the  vertical  shear  must  be  large  and  require  greater 
and  probably  deeper  compensating  flows  at  depth.   As  pre- 
viously discussed,  these  regions  will  most  likely  be  areas 
of  significant  deep  structure. 
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IV.   CONCLUSIONS 

A.   SUMMARY  OF  RESULTS 

The  distribution  of  temperature,  density,  heat,  and 
dynamic  heights  confirmed  that  the  phenomenon  of  temperature 
mesostructure  near  Arctic  ice  margins  was  related  to  the 
presence  of  an  ice  margin  and  the  associated  melting  pro- 
cesses.  Mesostructure  was  characterized  by  large  tempera- 
ture gradients  occurring  simultaneously  with  low  density 
gradients.   It  was  present  in  a  relatively  broad  zone, 
18  to  65  km  wide,  across  the  margins.   It  varied  consider- 
ably over  short  distances,  with  individual  elements 
generally  being  traceable  for  only  a  few  kilometers. 

The  assumption  of  a  generally  northward  flow  for  the 
area  of  study  was  supported  by  the  heat  excess  and  ice 
melt  content  calculations.   A  decrease  in  heat  content  and 

- 

a  corresponding  increase  in  the  ice  melt  content  of  the 
water  column  was  observed  as  one  crossed  the  ice  margin. 
These  changes  agreed  reasonably  well  enough  with  the  heat 
required  and  the  melt  water  that  would  have  been  formed  by 
the  melting  of  the  existing  ice  cover  to  support  the 
assumption  of  a  generally  northward  flow.   The  flow, 
however,  displayed  the  possibility  of  not  always  flowing 
directly  across  the  margins.   Crossings  2  and  6  displayed 
both  a  heat  and  an  ice  melt  content  decrease  as  one 
progressed  farther  under  the  ice.   This  could  have 
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indicated  a  flow  approaching  the  margin  obliquely  and  thus 
carrying  some  of  the  melt  water  away.   The  many  irregulari- 
ties also  present  in  both  the  heat  and  ice  melt  results 
suggest  the  presence  of  complex  flows  near  the  margins. 

Large  heat  exchanges  and  downward  mixing  throughout 
the  entire  water  column  were  evident  near  the  ice  margins. 
The  warm  surface  flow  was  cooled  at  the  surface  by  melting 
to  form  the  nose  feature.   The  cold,  dense  bottom  water, 
not  present  south  of  the  margin,  was  dissipated  by  downward 
mixing  in  the  water  column  creating  a  high  density  wedge- 
like feature.   The  broad  zones  of  low  vertical  density 
gradient  resulting  from  this  mixing  were  found  only  near 
the  margin  and  not  far  under  the  ice.   Thus,  the  creation 
of  mesostructure  by  the  turbulent  mixing  by  ice  keels 
as  suggested  by  Corse  (1974)  is  probably  a  doubtful 
mechanism  due  to  the  lack  of  evidence  of  mixing  farther 
under  the  ice. 

The  area  near  the  ice  margin  was  characterized  by  large 
temperature  gradients  along  isopycnals  and  by  the  irregularity 
and  tilting  of  the  isopycnals.   These  irregularities  and 
tilting  must  result  in  lateral  pressure  gradients  that 
evidently  are  balanced  by  complex  dynamic  flows  in  the 
marginal  area. 

Further  investigation  led  to  the  discovery  of  a  dynamic 
hill  which  extended  about  30  km  across  the  margin.   The  hill 
was  due  to  the  melting  process  and  was  notable  to  a  depth 
of  about  10  m.   It  was  interesting  that  the  dynamic  hills 
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were  about  the  same  height  as  the  depression  computed  by 
Neshyba  (1974)  for  the  freezing  process.   If  his  numbers 
are  correct,  the  dynamic  hill  for  the  melting  case  should 
be  many  times  larger  since  melting  occurred  much  more 
rapidly  than  freezing  and  was  concentrated  over  a  much 
narrower  width.   However,  it  appears  that  fresh  water  was 
lost  from  the  area  by  circulation,  thus  greatly  diminishing 
the  freshening  effect. 

The  magnitude  of  the  dynamic  hill  was  found  to  be  of 
lesser  importance  than  the  gradient  of  the  hill  near  the 
ice  margin.   The  lateral  pressure  gradient  which  resulted 
from  the  horizontal  gradient  in  dynamic  height  was  assumed 
to  be  balanced  by  internal  frictional  stresses,  and  when 
approximate  calculations  of  the  effect  of  this  balance  were 
made,  a  velocity  shear  at  the  surface  resulted.   This  velocity 
shear  was  applied  to  a  simple  velocity  profile  which  resulted 
in  a  deepening  of  the  flow  axis  and  an  acceleration  of  the 
flow  at  depth.   This  result  provides  a  plausible  explanation 
for  the  source  of  the  downward  mixing  that  is  observed 
at  the  margins.   This  hypothesis  is  also  supported  by  the 
correlation  of  larger  horizontal  dynamic  gradients  at  the 
surface  with  a  greater  sloping  of  the  deep  isopycnals  and 
the  presence  of  more  severe  deep  mesostructure .   The 
greater  dynamic  forces  at  the  surface  should  result  in 
larger  velocity  shears  and  thus  stronger  effects  at  depth. 

Shallow  and  deep  mesostructure  are  believed  to  be 
related  to  different  degrees  of  mixing  action.   Shallow 
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structure  was  correlated  with  the  presence  of  the  nose 
feature  while  deep  structure  was  related  to  large  dynamic 
slopes.   Deep  structure  was  usually  within  the  density  zone 
of  26.0-26.5  sigma-t  units  which  was  found  below  20  m 
depth.   The  deep  structure  was  located  in  the  vertically 
broadened  areas  of  this  density  zone.   The  low  density 
gradient  of  this  broadened  area  implies  mixing  processes, 
which  one  might  expect  from  the  large  gradients  in  dynamic 
height.   This  area  of  low  vertical  density  gradient  also 
was  one  where  interleaving  might  easily  occur  due  to  the 
small  amounts  of  energy  required  to  lift  a  parcel  of 
water.   Deep  structure  appeared  to  be  more  severe  when  the 
nose  feature  was  not  well  defined.   Thus  it  would  appear 
that  the  higher  dynamic  slopes  at  the  surface  cause  more 
severe  velocity  shears.   The  resulting  stronger  mixing 
effects  can  destroy  the  nose  feature  and  thus  explain  the 
usual  absence  of  the  nose  when  deep  structure  is  strong. 

Shallow  mesostructure  was  related  to  the  nose  feature 
and  was  generally  found  above  the  26.0  sigma-t  contour  and 
thus  usually  shallower  than  20  m  depth.   It  was  observed 
at  the  boundaries  of  the  nose  feature  where  water  types  of 
similar  density  and  different  temperatures  were  available 
to  interleave.   Shallow  structure,  with  two  exceptions, 
was  found  to  be  more  severe  in  the  cases  in  whicli  deep 
structure  was  weak  or  absent. 
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B.   RECOMMENDATIONS  TOR  FUTURE  INVESTIGATIONS 

As  a  result  of  this  study  it  is  evident  that  additional 
data  are  necessary  in  future  investigations.   The  probable 
variations  in  local  flows  near  the  ice  margin  make  it 
apparent  that  current  vector  data  are  needed  for  the  stations 
taken  near  the  margins.   Data  also  need  to  be  collected  far- 
ther into  the  ice  and  more  exact  data  on  ice  concentration 
and  thickness  would  be  helpful.   Time  series  data  on  the  STD 
parameters  and  currents,  taken  simultaneously,  would  also 
be  useful  in  determining  the  time  variance  of  the  phenomena. 
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APPENDIX  A 
STD  DATA  FOR  STATIONS  ANALYZED 

The  data  for  each  station  of  the  ice  crossings  analyzed 
follows.   Four  stations  appear  on  each  page  and  the  respec- 
tive station  number  appears  in  the  lower  part  of  each  plot. 
The  scales  for  the  horizontal  axes  are  located  at  the  top 
of  each  column  and  the  symbols  are  as  follows : 

ST,  sigma-t 

SV,  sound  velocity  (m  sec   ) 

S,  salinity  (°/oo) 

T,  temperature  (°C) 
The  symbol  for  each  of  the  above  parameters  is  at  the  bottom 
of  each  trace  and  an  S,  T,  or  V  (salinity,  temperature,  or 
sound  velocity)  appears  at  the  top  to  help  distinguish  indi- 
vidual traces.   The  gap  usually  located  at  about  10  depth 
on  each  trace  is  due  to  the  requirement  for  two  lowerings 
of  the  STD. 
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